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Abstract 

We present a new measurement of the W mass using the W — ► ev data 
from the D0 forward detectors at the Fermilab Tevatron pp Collider. This is 
the first measurement of the W mass with electron candidates in the range 
1.5 <| rj |< 2.5. We present measurements of the W mass using the transverse 
mass, the electron transverse momentum and the neutrino transverse momen- 
tum, and the combined result using all three techniques. The combination of 
the forward detector measurement with the previous measurements using the 
central detector gives a new precise measurement of the W mass from D0. 
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In the standard model of the electroweak interactions (SM) [T|], the mass of the W boson 
is predicted to be 



M _ Ara(Mj) y 1 

\ y/2G F ) sin ^ VI- ' 

In the "on-shell" scheme @ cos# w = My//Mz, where Mz is the Z boson mass. A measure- 
ment of My/, together with Mz-, the Fermi constant {Gf), and the electromagnetic coupling 
constant (a), determines the electroweak radiative corrections Ary/ experimentally. Purely 
electromagnetic corrections are absorbed into the value of a by evaluating it at Q 2 = M§ [|3] . 
The dominant contributions to Ary/ arise from loop diagrams that involve the top quark 
and the Higgs boson. If additional particles which couple to the W boson exist, they give 
rise to additional contributions to Ary/- Therefore, a measurement of My/ is one of the 
most stringent experimental tests of SM predictions. Deviations from the predictions may 
indicate the existence of new physics. Within the SM, measurements of My/ and the mass 
of the top quark constrain the mass of the Higgs boson. 

This paper reports a new measurement of the W boson mass using electrons detected in 
forward calorimeters. We used data recorded by the D0 detector || during the 1994-1995 
run of the Fermilab Tevatron pp collider. This forward electron measurement complements 
our previous measurement with central electrons || because the more complete rapidity cov- 
erage gives useful constraints on model parameters that permit reduction of the systematic 
error, in addition to increasing the statistical precision. Combining this measurement with 
our previously published measurements with central electrons using data taken in 1992-1993 
and 1994-1995, determines the W boson mass to a precision of 93 MeV. 

At the Tevatron, W bosons are produced mainly through qq annihilation. We detect 
them by their decays into electron-neutrino pairs, characterized by an isolated electron || 
with large transverse momentum (p?) and significant transverse momentum imbalance {^t)- 
The p't is due to the neutrino which escapes detection. Many other particles of lower 
momenta, which recoil against the W boson, are produced in the breakup of the proton and 
antiproton. We refer to them collectively as the underlying event. 

At the trigger level we require j/y > 15 GeV and an energy cluster in the electromag- 
netic (EM) calorimeter with pt > 20 GeV. The cluster must be isolated and have a shape 
consistent with that of an electron shower. 

During event reconstruction, electrons are identified as energy clusters in the EM 
calorimeter which satisfy isolation and shower shape cuts and have a drift chamber track 
pointing towards the cluster centroid. We determine forward electron energies by adding the 
energy depositions in the first ~ 40 radiation lengths of the calorimeter in a cone of radius 
20 cm, centered on the highest energy deposit in the cluster. Fiducial cuts reject electron 
candidates near calorimeter module edges and ensure a uniform calorimeter response for the 
selected electrons. The electron momentum (p(e)) is determined by combining its energy 
with its direction which is obtained from the shower centroid position and the drift chamber 
track. The trajectories of the electron and the proton beam define the position of the event 
vertex. 

We measure the sum of the transverse momenta of all the particles recoiling against the 
W boson, u T = Ei sin 6^, where Ei is the energy deposition in the i th calorimeter cell 
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and 6{ is the angle denned by the cell center, the event vertex, and the proton beam. The 
unit vector Ui points perpendicularly from the beam to the cell center. The calculation of ut 
excludes the cells occupied by the electron. The sum of the momentum components along 
the beam is not well measured because of particles escaping through the beam pipe. From 
momentum conservation we infer the transverse neutrino momentum, priy) = — Pr(e) — Wr, 
and the transverse momentum of the W boson, pt{W) = —ut- 

We select a W boson sample of 11,090 events by requiring priy) > 30 GeV, ut < 15 GeV, 
and an electron candidate with 1.5 < \r]\ < 2.5 and pr(e) > 30 GeV. 

Since we do not measure the longitudinal momentum components of the neutrinos 
from W boson decays, we cannot reconstruct the ev invariant mass. Instead, we ex- 
tract the W boson mass from the spectra of the electron pt and the transverse mass, 
rriT = ^/2pT(e)pT(i / )(l — cos A0), where A0 is the azimuthal separation between the two 
leptons. We perform a maximum likelihood fit to the data using probability density func- 
tions from a Monte Carlo program. Since neither my nor pr(e) are Lorentz invariants, we 
have to model the production dynamics of W bosons to correctly predict the spectra. The 
rriT spectrum is insensitive to transverse boosts at leading order in px{W) / Myy and is there- 
fore less sensitive to the W boson production model than the Pt{&) spectrum. On the other 
hand, the rriT spectrum depends strongly on the detector response to the underlying event 
and is therefore more sensitive to detector effects than the pr(e) spectrum. The shape of 
the neutrino priy) spectrum is sensitive to both the W boson production dynamics and the 
recoil momentum measurement. By performing the measurement using all three spectra we 
provide a powerful cross check with complementary systematics. 

Z bosons decaying to electrons provide an important control sample. We use them to 
calibrate the detector response to the underlying event and to the electrons, and to constrain 
the model for intermediate vector boson production used in the Monte Carlo simulations. 

A Z — > ee event is characterized by two isolated high-pr electrons. We trigger on events 
with at least two EM clusters with p T > 20 GeV. We accept Z — ► ee decays with at least 
one forward electron in the pseudorapidity range 1.5 < \rj\ < 2.5, where rj — — In (tan|), 
and the other electron to be either forward or central with pseudorapidity \rj\ < 1.0. The 
central electron is required to have px > 25 GeV but is allowed not to have a matching 
drift chamber track. The forward electron candidate is required to have px > 30 GeV and 
a matching drift chamber track. This selection accepts 1,687 events. 

For this measurement we used a fast Monte Carlo program developed for our central 
electron analysis ||, with some modifications in the simulation of forward electron events. 
The program generates W and Z bosons with the rapidity and px spectra given by a calcu- 
lation using soft gluon resummation [0] and the MRSA |T3| parton distribution functions. 
The line shape is a relativistic Breit-Wigner, skewed by the mass dependence of the par- 
ton luminosity. The measured intrinsic widths [p|,[T0|] are used. The angular distribution of 
the decay electrons includes a pr(W)-dependent 0(a s ) correction fCTf. The program also 



generates W — > ez/7 [0, Z — > eej |[L2|1 , and W — > rv — >• evvv decays. 

The program smears the generated p(e) and ut vectors using a parameterized detec- 
tor response model and applies inefficiencies introduced by the trigger and event selection 
requirements. The model parameters are adjusted to match the data and are discussed 
below. 

The electron energy resolution is described by sampling, noise, and constant terms. In 



6 



the Monte Carlo simulation of forward electrons we use a sampling term of 15.7%/ ^/p/GeV, 
derived from beam tests. The noise term is determined by pedestal distributions derived 
from the W data sample. We constrain the constant term to cem = 1-0^?;q% by requiring 
that the width of the dielectron invariant mass spectrum predicted by the Monte Carlo 
simulation is consistent with the Z data (Fig. [I]). 




70 75 80 85 90 95 100 105 110 

m(ee)(GeV) 

FIG. 1. The dielectron invariant mass distribution of the Z data (•). The solid line shows 
the fitted signal plus background shape and the small shaded area the background. The arrows 
indicate the fit window. 



Beam tests show that the electron energy response of the calorimeter can be parameter- 
ized by a scale factor «em and an offset 5em- We determine these in situ using Z — > ee decays. 
We obtain for forward electrons <5em = —0.1 ± 0.7 GeV and «em = 0.95179 ± 0.000187 by 
fitting the observed mass spectra while constraining the resonance masses to the measured 
value of the Z boson mass [ |Tu||i"5| . The uncertainty in «em is dominated by the finite size 
of the Z sample. Figure p] shows the observed dielectron mass spectrum from the dielectron 
sample, and the line shape predicted by the Monte Carlo simulation for the fitted values of 
cem, «emj an d <5em- 

We calibrate the response of the detector to the underlying event, relative to the EM 
response, using Z boson data sample. In Z — > ee decays, momentum conservation requires 
Pr(ee) = —Ut, where pr(ee) is the sum of the two electron pt vectors. To minimize sensi- 
tivity to the electron energy resolution, we project ut and pr(ee) on the inner bisector of 
the two electron directions, called the ry-axis (Fig. §). We call the projections u v and p v {ee). 



ft 






U, j 






p,(ee) 1 












FIG. 2. The definition of the 77-axis (left). The plot of u^/Rr 
(•) and simulation ( — ). 
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+ p v (ee) (right) for the data 



Detector simulations based on the geant program WM predict a detector response to 
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the recoil particle momentum of the form R rec = a Tec + /3 rec \n(pT/GeV). We constrain 
a rec and /3 rec by comparing the mean value of u v + p v (ee) with Monte Carlo predictions for 
different values of the parameters. We use a rec = 0.693 ± 0.060 and (3 rec = 0.040 ± 0.021 
with a correlation coefficient of —0.98 obtained in central electron analysis ||. We check 
that Z events with both electrons in the end calorimeter give a recoil response measurement 
statistically consistent with the above (Figure §). 

The recoil momentum resolution has two components. We smear the magnitude of 
the recoil momentum with a resolution of s rec / ^Jpt/GqV . We describe the detector noise 
and pile-up, which are independent of the boson p? and azimuthally symmetric, by adding 
the $t from a random pp interaction, scaled by a factor a m b, to the smeared boson p T . 
To model the luminosity dependence of this resolution component correctly, the sample of 
pp interactions was chosen to have the same luminosity spectrum as the W sample. We 
constrain the parameters by comparing the observed rms of li^/Rrec + P»?(ee) with Monte 
Carlo predictions. We use s rec = 0.49 ± 0.14 and a m b = 1.032 ± 0.028 with a correlation 
coefficient of —0.60 measured in central electron analysis |J. Since we exlude the cells 
occupied by the electrons, the average transverse energy flow, St = ^i?jsin#j, is higher 
for the W sample than for the Z sample. This bias is caused by requiring the identification of 
two electrons in the Z sample versus one in the W sample. The larger energy flow translates 
into a slightly broader recoil momentum resolution in the W sample. a m b is corrected by a 
factor 1.03 ± 0.01 to account for this effect in the W boson model. The p n balance width is 
in good agreement between data and Monte Carlo for our Z event sample. Hence we use the 
same recoil resolution for forward electorn W events as for the central W events || . Figure |2| 
shows a plot of Mjj/Rrec + p v {ee) when both electrons are in the forward calorimeters. 

Backgrounds in the W sample are W — > rv — > evvv decays (1.0%, which is included into 
the Monte Carlo simulation ), hadrons misidentified as electrons (3.64%±0.78%, determined 
from the data ), and Z -> ee decays (0.26%±0.02%, determined from HERWIG/GEANT 
simulations). Their shapes are included in the probability density functions used in the fits. 




FIG. 3. Spectra of (a) mj, (b) pr{e) and (c) pr{v) from the data (•), the fit ( — ), and the 
backgrounds (shaded). The arrows indicate the fit windows. 

In principle, if the acceptance for the W — > ev decays were complete, the transverse mass 
distribution or the lepton p T distribution would be independent of the W rapidity. However, 
cuts on the electron angle in the laboratory frame cause the observed distributions of the 
transverse momenta to depend on the W rapidity. Hence a constraint on the W rapidity 
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distribution is useful in contraining the production model uncertainty on the W mass. We 
introduce in the Monte Carlo a scale factor k v for the W boson rapidity as r] W — > k v r]w 
and extract it by fitting Monte Carlo electron rapidity distribution simulated with different 
value of k v to the data. The scale factor is found to be consistant, within errors, with unity 
with uncertainty 1.6%. 

The fit to the txit distribution (Fig. 0(a)) yields Mw = 80.766 GeV with a statistical 
uncertainty of 108 MeV. A x 2 test gives \ 2 = 17 for 25 bins which corresponds to a confidence 
level of 81%. The fit to the p T (e) distribution (Fig. 0(b)) yields M w = 80.587 GeV with 
a statistical uncertainty of 125 MeV. The confidence level of the y 2 test is 5%. The fit to 
the Pt(v) distribution (Fig. 0(c)) yields M w = 80.726 GeV with a statistical uncertainty of 
163 MeV and confidence level 33%. 

We estimate systematic uncertainties in Mw from the Monte Carlo parameters by varying 
them within their uncertainties (Table I). In addition to the parameters described above, the 
calibration of the electron polar angle measurement contributes a significant uncertainty. We 
use muons from pp collisions and cosmic rays to calibrate the drift chamber measurements 
and Z — > ee decays to align the calorimeter with the drift chambers. Smaller uncertainties 
are due to the removal of the cells occupied by the electron from the computation of Ht, the 
uniformity of the calorimeter response, and the modeling of trigger and selection biases. 



TABLE I. Uncertainties in the W boson mass measurement in MeV. 



Source 


Forward 


Forward + Central 


W boson statistics 


107 


61 


Z boson statistics 


181 


59 


Calorimeter linearity 


52 


25 


Calorimeter uniformity 8 


Electron resolution 


42 


19 


Electron angle calibration 


20 


10 


Recoil response 


17 


25 


Recoil resolution 


42 


25 


Electron removal 


2 


12 


Selection bias 


5 


3 


Backgrounds 


20 


9 


Parton distribution functions 


35 


15 


Parton luminosity 


2 


4 


Pt(W) spectrum 


25 


15 


W boson width 


10 


10 


radiative decays 


1 


12 



The uncertainty due to the model for W boson production and decay consists of several 
components (Table I). We assign an uncertainty that characterizes the range of variations 
in M\y obtained when employing several recent parton distribution functions: MRST ||, 
MRSA jnj, MRSR2 [T£j, CTEQ5M |T7|, and CTEQ4M |T8J. We allow the p T (W) spec- 
trum to vary within constraints derived from the pr(ee) spectrum of the Z data and from 
TBI 



A 



QCD 



The uncertainty due to radiative decays contains an estimate of the effect of 
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neglecting double photon emission in the Monte Carlo simulation [19 



The fit to the my spectrum results in a W boson mass of 80.766 ± 0.108(stat) ± 
0.208(syst) GeV, the fit to the p T (e) spectrum results in 80.587 ± 0.125(stat) ± 
0.217(syst) GeV, and the fit to the priy) spectrum results in 80.726 ± 0.163(stat) ± 
0.307(syst) GeV. The good agreement of the three fits shows that our simulation mod- 
els the W boson production dynamics and the detector response well. Fits to the data in 
bins of luminosity, 0(e), r](e), and ut do not show evidence for any systematic biases. 

We combine all the six measurements ( fits to central electron || and forward electron 
W boson events using three techniques ). We obtain the combined 1994-1995 measurement 
M w = 80.487 ± 0.096 GeV. The x 2 is 4.6/5 dof, with a probability of 46%. Combining with 
the measurement from the 1992-93 data gives the 1992-95 data measurement of M\y = 
80.474 ± 0.093 GeV. 
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FIG. 4. Direct W boson and top quark mass measurements by the D0 pO] experiment. The 



bands show SM predictions for the indicated Higgs masses [21j 



Using Eq. [T] we find Ar\y = —0.0317 ± 0.0061, which establishes the existence of elec- 
troweak corrections to M\y at the level of five standard deviations. Figure ^ compares the 
direct measurements of the W boson and top quark masses to SM predictions. 
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and acknowledge support from the Department of Energy and National Science Foundation 
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of Atomic Energy and Science and Education (India), Colciencias (Colombia), CONACyT 
(Mexico), Ministry of Education and KOSEF (Korea), and CONICET and UBACyT (Ar- 
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